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DOES THE SLIM-DISK MODEL CORRECTLY CONSIDER PHOTON-TRAPPING EFFECTS? 
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ABSTRACT 

We investigate the photon-trapping effects in the super-critical black hole accretion flows by solving 
radiation transfer as well as the energy equations of radiation and gas. It is found that the slim-disk 
model generally overestimates the luminosity of the disk at around the Eddington luminosity (Le) and is 
not accurate in describing the effective temperature profile, since it neglects time delay between energy 
generation at deeper inside the disk and energy release at the surface. Especially, the photon-trapping 
effects are appreciable even below L ~ Le, while they appear above ~ 3Le according to the slim disk. 
Through the photon-trapping effects, the luminosity is reduced and the effective temperature profile 
becomes flatter than r -3 / 4 as in the standard disk. In the case that the viscous heating is effective only 
around the equatorial plane, the luminosity is kept around the Eddington luminosity even at very large 
mass accretion rate, M 3> Le/c 2 . The effective temperature profile is almost flat, and the maximum 
temperature decreases in accordance with rise in the mass accretion rate. Thus, the most luminous radius 
shifts to the outer region when Mj (Le/c 2 ) » 10 2 . In the case that the energy is dissipated equally at any 
heights, the resultant luminosity is somewhat larger than in the former case, but the energy-conversion 
efficiency still decreases with increase of the mass accretion rate, as well. The most luminous radius stays 
around the inner edge of the disk in the latter case. Hence, the effective temperature profile is sensitive 
to the vertical distribution of energy production rates, so is the spectral shape. Future observations of 
high L/Le objects will be able to test our model. 

Subject headings: accretion: accretion disks — black hole physics — radiative transfer 



1. INTRODUCTION 

X-ray binary sources and active galactic nuclei (AGNs) 
emit enormous energy in radiation, and it is believed that 
the accretion disks are the main place of energy release 
in those objects. The standard-disk model was proposed 
by Shakura & Sunyaev (1973) as a very efficient con- 
version mechanism of gravitational energy of accreting 
gas into radiation energy (Pringle 1981; Frank, King, & 
Raine 1985; Kato, Fukue, & Mineshige 1998). The stan- 
dard model has been very successful in describing opti- 
cally thick flow structure, as long as the mass accretion 
rate, M, is less than the critical mass accretion rate, 
M CT it = Le/c 2 , with Le being the Eddington luminos- 
ity given by At: cG M m v / a^, where c is the light velocity, 
M is the black hole mass, m p is the proton mass, and <7t 
is the Thomson scattering cross-section. However, sub- 
critical accretion is not always guaranteed, since the mass 
accretion rate is determined by something other than the 
central star itself, i.e., companion star in binary system 
(Alme k Wilson 1976; Spruit & Ritter 1983; King et al. 
1997; Koyama et al. 1999) or circumnuclear star clusters 
in AGNs (Norman & Scoville 1988; Umemura, Fukue, & 
Mineshige 1998; Ohsuga et al. 1999). Therefore, it is no 
wonder that the mass accretion rate can greatly exceed 
the critical value. 

When the mass accretion rate is comparable to or more 



than the critical limit, the disk becomes radiation pres- 
sure dominant. In the case of M>10M cr it, moreover, the 
disk is moderately geometrically thick and advective en- 
ergy transport becomes substantial. Such a disk is called a 
slim disk and has been investigated in detail (Abramowicz 
et al. 1988; Szuszkiewicz, Malkan, & Abramowicz 1996; 
Wang et al. 1999; Watarai & Fukue 1999; Watarai et al. 
2000). The unique feature in high M/M cr j t systems is that 
photons are trapped (Katz 1977; Begelman 1978; Begel- 
man & Meier 1982 ; Flammang 1984; Blondin 1986; Colpi 
1988; Wang & Zhou 1999 ). In optically thick accretion 
flow, frequent interaction between matter and photons de- 
lays liberation of the radiation energy arising at deeper 
inside the disk. Thus, the radiation energy is trapped in 
the flow and advected inward. Such a photon-trapping 
plays an important role when the radiative diffusion time- 
scale is longer than the accretion time-scale. Since the 
trapped radiation energy can fall onto the black hole with 
the accreting gas without being radiated away, the ob- 
served luminosity is reduced in the black-hole accretion 
flow. In contrast, the advected energy should be finally 
radiated at the stellar surface in the case of neutron star 
(Houck & Chevalier 1992). 

While most studies focused on spherical accretion case 
with some discussion on Comptonization, we here focus 
on the disk accretion case. In the slim-disk approach, the 
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vertically-integrated approximation is adopted; namely 
the radiative flux at the disk surface is related to the tem- 
perature on the equatorial plane, T c , through the relation, 
F s z urf <~ crT c 4 /r, where a is the Stcfan-Boltzmann constant 
and t is the vertical Thomson scattering optical depth for 
a half of the disk. However, this holds only when the radia- 
tive diffusion time-scale is shorter than the accretion time- 
scale. In other words, the photon-trapping effects are not 
exactly taken into consideration in the slim-disk formula- 
tion. Eggum, Coroniti, & Katz (1987, 1988) initiated the 
two-dimensional radiation-hydrodynamical simulations by 
assuming the equilibrium between gas and radiation. Re- 
cently, the improved simulations, in which the energy of 
gas and radiation are separately treated, were performed 
by Okuda, Fujita, & Sakashita (1997), Fujita & Okuda 
(1998) and Kley & Lin (1999). 

In this work, by solving radiation transfer, we investi- 
gate the photon-trapping effects, paying attention to en- 
ergy transmission inside the disk in the optically-thick 
black hole accretion flows. The validity of the slim-disk 
model is also discussed in terms of the luminosity and 
the effective temperature profile. In §2, we present ba- 
sic considerations of the photon-trapping effects based on 
the comparison of the accretion time-scale and the radia- 
tive diffusion time-scale. Model and basic equations for 
numerical simulations are given in §3, and we describe re- 
sults in §4. Finally, §5 and §6 are devoted to discussion 
and conclusions. 

2. BASIC CONSIDERATIONS 

The radiation energy generated near the equatorial 
plane diffuses toward the disk surface at the speed of 
~ c/3r (Mihalas & Mihalas 1984), so that the time-scale 
of radiative diffusion is tdiff = H/(c/3t), where H is the 
disk half-thickness. Since the accretion time-scale, i acc , is 
given by —r/v r , the condition that the radiation energy in 
the disk is trapped in the flow and falls onto the black hole 
is written as H/(c/3t)> — r/v r , where r is the radius and 
v r is the accretion velocity of the flow. Using the relation, 
r = (7TS/2mp, and the continuity equation of the accre- 
tion disk, M — — 27rru r £, we obtain the photon-trapping 
radius as 

3 

rtrap = ^rhr s h, (1) 

where £ is the surface density of the disk, m is the mass ac- 
cretion rate normalized by the critical mass accretion rate, 
fa = M I Merit, f g is the Schwarzschild radius, and ft is the 
ratio of the half disk-thickness to the radius ft = H/r, 
respectively. Hence, the photon-trapping effects perform 
significant role in the super-critical accretion flow when 
»Vap > 3r g , namely at fn>2, since ft is of the order of 
unity in radiation pressure dominant region. 

The energy dissipated by viscosity at the regions of 
f < ^trap can not be released but is carried onto the black 
hole, as long as the radiative diffusion is the predominant 
process for energy transport. Since the radiation arising 
merely at r > r trap can pass through the disk body to go 
out from the surface, the luminosity is roughly estimated 

by 
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where r; n is the inner edge of the disk, Q v j s is the half 
vertically-integrated viscous heating rate, and the com- 
pressional heating is neglected. Adopting the viscous hat- 
ing rate as 
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(Shakura & Sunyaev 1973; Lynden-Bell & Pringle 1974), 
with fl K being the Keplerian angular speed, we have in 
Newtonian approximation 



for m <C 2n n /(3r g h), and 



L ~ -M cxKt c z 



2 

3ft 



4n, 

2 
3 



2 

3ft 



3/2 



r„m 



1/2 



(4) 



(5) 



if fa is larger than 2ri n /(3r g ft). It is found in the limit of 
large fa (s> 1) that the luminosity is roughly constant ir- 
respective of m so that the energy-conversion efficiency, 
■n = L/Mc 2 ~ M crit /2Mft = (2ftm)-\ should be re- 
markably reduced in the super-critical accretion disk with 
m 3> 1. Also, it is important to note that the photon- 
trapping effects do not depend on the accretion velocity 
nor viscosity for a fixed m, since the accretion velocity in- 
creases if a rises, but at the same time, the radiative dif- 
fusion velocity also increases in accordance with decrease 
in the surface density of the disk (i.e., S oc v^ 1 oc a -1 ). 

So far, we have considered the limiting model that the 
viscous heating occurs only in the vicinity of the equato- 
rial plane. Next, we investigate the photon-trapping ef- 
fects for the model that the gas is heated up in proportion 
to the density. Just as in the previous model, the photon- 
trapping effects are appreciable within the trapping radius, 
however, the radiation arising near the disk surface is not 
trapped in the flow. The condition for the photon trapping 
to occur at z is given by {H — z) / (c/ 3t z ) > —r / v r ■, where z is 
the vertical height and t z is the vertical optical depth mea- 
sured from the surface of the disk, t z = J (p<jt /m v ) dz, 
with p being the gas density. Assuming that the gas 
density is constant in the vertical direction for simplic- 
ity, namely t z = t(H — z)/H, and using the continuity 
equation of the accretion disk as well as equation (1), we 
find that the energy dissipated at t z < (r/rt rap ) 1 / 2 T can be 
radiated away even within the trapping radius. Therefore, 
the luminosity is estimated by 
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, (6) 



when fa > 2r- m /{3r g h). In this model, the luminosity is 
not constant but depends on the mass accretion rate as 
oc m 1 / 2 at m > 1. The energy-conversion efficiency de- 
creases as J] oc m -1 / 2 . If fa < 2n n /(3r g ft), the photon- 
trapping effects do not appear and the observed luminos- 
ity is given by equation (4). Here, it is stressed that the 
assessment with using the radiative diffusion velocity is 
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valid when the most of radiation energy arises at deeper 
inside the disk. If the viscous heating occurs only around 
the disk surface, the photon trapping would be ineffective. 
The reason of this is that the photons arising at near disk 
surface (r 2 <C 1) can easily escape from the flow since an 
average number of the scatterings is ~ t z , in spite of r z for 
photons emitted at t z 3> 1 (Sunyaev & Titarchuk 1985). 

In Figure 1, we show the expected luminosity changes 
as functions of to assessed from equations (4), (5), and (6). 
To compare with the numerical results (§4.2), the former 
and latter models are represented as the analytical models 
A (thick solid curve) and B (dashed curve), respectively, 
and h = 0.5 as well £lS f^iYl = 3rv are assumed. In the 
sub-critical regime, the luminosities for both models vary 
along the thin solid line. We also display the luminosity of 
the slim disk (dotted curve). These two models are both 
extreme in the opposite sense, and actual situation will 
lie just in between. We thus expect that the luminosity 
increases moderately but rises slowly than to 1 / 2 . Hence, 
the slim disk is expected to overestimate the luminosity 
at L>Le, namely m> a few. In §4.2, we discuss about 
behavior of the luminosity in more detail. To assess the 
radiative flux as a function of radius, we need to calculate 
the radius, at which the energy emerged at outer region is 
released. Note that finite disk size effect enhances photon 
trapping, and the disks would be even fainter (see §5.2). 
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Fig. 1. — The luminosity is plotted against the mass accretion 
rate, where the luminosity and mass accretion rate are normalized 
by the Eddington luminosity and the critical mass accretion rate, 
respectively. The thick solid and dashed curves are the luminosi- 
ties analytically estimated in §2. The circles and squares are the 
numerical results for models A and B, respectively. The luminosity 
should increase along the thin solid line, if energy-conversion effi- 
ciency is constant. It is found, however, that the energy-conversion 
efficient is on decrease in the super-critical accretion flow. Clearly, 
the slim disk (dotted curve) considerably overestimates the luminos- 
ity as compared with model A at whole super-critical regime and 
with model B when m ~ 10 — 10 3 . The numerical results are roughly 
consistent with analytical ones. 

3. MODEL AND BASIC EQUATIONS FOR NUMERICAL 
SIMULATIONS 



In the previous section, we roughly estimate the photon- 
trapping effects by comparing the accretion time-scale and 
the radiative diffusion time-scale. In this section, we ex- 
amine the photon-trapping effects by practically solving 
radiation transfer, and the energy equations of gas and 
radiation. Since we are more concerned with the photon- 
trapping effects themselves rather than flow structure, we 
employ a simple model for the accretion flow. Here, the 
cylindrical coordinate, (r, tp, z), is used. We consider that 
the accretion disk is axisymmetric and steady in the Eule- 
rian description; d/dip = d/dt = 0. The radial component 
of the velocity is expressed in terms of free-fall velocity as 

'GM\ 1/2 



(7) 



where £ is a constant parameter, and its vertical compo- 
nent is prescribed as 

V z = -V r , (8) 

r 

i.e., we assume convergence flow. Note that v r is related to 
the viscosity parameter a through v r ~ a (H/r) c s , where 
c s is the sound velocity. We focus on the accreting ring 
element whose geometrical width and thickness are Ar 
and 2H, respectively, where Ar <C r (see Figure 2). We 
suppose the structure of the accretion disk to be locally 
plane parallel, that is to say, each ring is composed of N 
layers in the z-direction, and the thickness of each layer is 
Az = 2H/N. We solve the time-dependent energy fields of 
gas and radiation in the ring element during the course of 
accretion motion until the element reaches the inner edge 
of the disk. Wc adopt uniform density profile, 
E(r) 



p(r, z) 



const, m z, 



(9) 



2H(r) 

and consider both of the viscous heating and compressional 
heating due to converging inflow. Since the flow is steady 
(d/dt — 0), time coordinate can be transformed to the spa- 
tial coordinates; i.e., D/Dt = v r d/dr+v z d/dz with v r and 
v z being given by equations (7) and (8). We thus express 
the radiative flux at the disk surface as a function of radius 
and then calculate the resultant luminosity of the disk. 



ring elements 




black hole accretion disk 



Fig. 2. — Schematic view of the black-hole accretion system ex- 
plaining our calculation method. We calculate the evolution of the 
energy fields in a moving ring clement that consists of N layers, un- 
til the element reaches the inner edge of the disk, assuming layers 
in each ring to be locally plane-parallel. Since the steady flow in 
the Eulerian description is considered, we set d/dt = 0, and thus 
replace D/Dt = v r 8/dr + v z d/dz. 

In plane parallel approximation, the radial and az- 
imuthal components of the radiative flux arc null. Also, 
non-diagonal components of the radiation stress tensor are 
null. Hence, using equations (7) and (8), we write the en- 
ergy equations of radiation and gas as 
d z d \ ( 'E s 
r dz 
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dF z v r 
dz r 



+4ttkB - ckE, 



(10) 



and 



or r dz 



3 v r 

--—p gas -4irKB+CK,E+q v 
2 r 



(11) 

respectively, where E is the radiation energy density, F z 
is the radiative flux in the z-direction, P" is the diagonal 
components of the radiation stress tensor, k is the absorp- 
tion coefficient, B is the blackbody intensity written as 
B = (rTgas/TT with T gas being the gas temperature, e is 
the internal energy density, p gas is the gas pressure, and 
Qvis is the viscous heating rate per unit volume (Mihalas & 
Klein 1982; Mihalas & Mihalas 1984; Fukue, Kato, & Mat- 
sumoto 1985; Stone, Mihalas, & Norman 1992). To close 
the set of equations, we apply the flux-limited diffusion 
approximation for the radiative flux and stress tensor, 

CX9E (12) 



F* = - 



X dz ' 
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p<P<P = _(!_/)£, 



(13) 



and 



P zz = fE, (14) 

with x being the extinction coefficient, where A, R, and / 
are defined as 

A - ^4±^, (15) 
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/ = A + A 2 i? 2 (17) 

(Turner & Stone 2001). This approximation holds both in 
the optically thick and thin regimes. In the optically thick 
limit, we find A — > 1/3 and / — > 1/3 because of R — ► 0. 
In the optically thin limit of R — > oo, on the other hand, 
we have \F Z \ = cE, P rr = P^ = 0, and P zz = E. These 
give correct relations in the optically thick diffusion limit 
and optically thin streaming limit, respectively. Thus, we 
can rewrite the radiation energy equation (10) as 
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These nonlinear equations (11) and (18) are integrated 
iteratively by the Ncwton-Raphson method with the 
Gauss- Jordan elimination for a matrix inversion, by cou- 
pling with the continuity equation, 



M = -27TTO r S, 



(19) 



the density profile given by equation (9), and the equation 
of state, 

2 



Pgas 
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As to the vertical distribution of viscous heating rates, we 
consider two extreme models. In model A, we assume that 



the viscous heating is effective only in the vicinity of the 
equatorial plane at \z\ < 10~ 2 H, that is to say, 



Qvisiz) 



_ / l0 2 Qvis/H \z\ < 1Q- 2 H 
\z\ > 10~ 2 H 



(21) 



where Q v is is given by equation (3). In model B, we as- 
sume that the gas is heated up uniformly, independently 
of z, that is, 



/ Qvis 

qvis(z) = —j-j- = const. 



(22) 



Throughout the present study, the black hole mass is fixed 
to be 10M Q and the inner edge of the disk is taken to be 
at ri n = 3r g . Moreover, we assume H to be 0.5r (i.e., 
h = 0.5), and N — 100 is employed. This assumption of 
H = 0.5r would be reasonable in the super-critical accre- 
tion regime, since H/r of the radiation pressure dominant 
disk becomes of the order of unity. In the sub-critical 
regime, on the other hand, H would be much smaller than 
r, but then the most of emergent energy can be released 
immediately without being trapped. Thus, our numerical 
simulations are valid in both super-critical and sub-critical 
accretion regimes for evaluating the radiative flux at the 
disk surface. 



4. RESULTS 

4.1. Reduction of the radiative flux 

The radiative flux is reduced due to the photon trap- 
ping in the vicinity of the inner edge of the disk when 
m 3> 1. Figure 3 shows the ratios of the radia- 
tive flux at the disk surface, -F„ r f, to the vertically- 
integrated viscous heating rate for several m as func- 
tions of radius. The solid and dashed curves correspond 
to model A and model B, respectively. Here, we con- 
sider the situation that the Thomson scattering is pre- 
dominant over the absorption, n = 10~ 2 pcr^/m p . We 
set £ in equation (7) for different m to be (m, £) = 

(0.1, io- 4 ), (1, io- 3 ), (10, io- 3 ), (10 2 , io- 2 ), (10 3 , io- 2 ), 

and (10 4 ,0.1) in order for the disk to be optically thick 
for Thomson scattering. Then, the viscosity parameter a 
is of the order of 10~ 4 in £ = 10~ 4 and 0.1 in £ = 0.1, 
since a ~ £//i 2 from equations (7) and (8). As shown in 
Figure 3, the ratio of F z mf /Q v i s is around unity when the 
accretion rate is sub-critical, rh<l. In the super-critical 
accretion disk, rh » 1, on the other hand, the ratio is 
much smaller than unity in the vicinity of the inner edge 
of the disk, but is still close to unity at the outer region. 
The ratio rises at r ~ 3r 5 again, but this rise is caused 
by the decline of Q v ; s near the inner boundary. The ab- 
solute value of the radiative flux is still suppressed there. 
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Fig. 3. — The ratios of radiative flux at the disk surface to the 
vertically integrated viscous heating rate as functions of radius for 
rh = 0.1, 10, 10 2 , 10 3 , and 10 4 . Here, the radius is normalized by the 
Schwarzschild radius. The solid and dashed curves indicate models 
A and B, respectively. By the photon-trapping effects, the ratios 
are much smaller than unity at the regions of r < rtrap in the case 
of m > 1, where rtrap is the photon-trapping radius given by equa- 
tion (1). It is also found that the photon-trapping effects arc more 
conspicuous in model A than in model B. If rh < 1, on the other 
hand, we found ^g UT f /Qvia ~ 1, as in the standard disk. 

This reduction of the radiative flux is due to the pho- 
ton trapping, and it can be understood by the time delay 
owing to energy transport in the optically thick medium. 
The emergent radiation energy is diffused at the speed 
of ~ c/3t, and arrives at the disk surface after t^is- If 
this time delay is comparable to or longer than the accre- 
tion time-scale, idiff >£acc, the generated radiation energy 
is advected inward, so that the ratio of F s z urf /Q v is should 
deviate from unity. To conclude, the energy dissipated by 
viscosity can be released from the surface immediately at 
the same radius, as in the standard disk, only in the outer 
region for rh>l but in the whole region for rh < 1, where 
the radiative diffusion time-scale is much shorter than the 
accretion time-scale. 

As shown in Figure 3, the photon-trapping effects are 
appreciable within mr g in both models A and B. This 
is good agreement with the analytical prediction of r trap , 
[see equation (1)]. However, the photon-trapping effects 
are more conspicuous in model A than in model B, since 
in model B the energy dissipated near the disk surface 
can still escape from the accretion flow even within the 
trapping radius. In model A, in contrast, most of energy 
is dissipated in the vicinity of the equatorial plane and 
is trapped in the flow. Realistic situation will be just in 
between. 

Here, it should be stressed again that the photon- 
trapping effects themselves depend on rh and are indepen- 
dent of v r when rh is fixed, as was mentioned in §2, unless 
the disk is optically thin. To check if this is the case, we 
examine the cases with different £ values; £ = 1CP 3 and 
0.1, besides £ = 10 -2 , for rh = 10 2 , finding no changes. 



[However, if the absorption coefficient is very small, the 
luminosity would decrease in the case of rapid accretion, 
since thermal energy of the gas can not be converted into 
radiation energy (see §4.6).] 

4.2. Luminosity 

The luminosity of the disk is significantly reduced by 
the photon-trapping effects. In Figure 1, the calculated 
luminosity is also plotted against rh with symbols. Here, 
the circles and squares indicate the luminosities for mod- 
els A and B, respectively. The thin solid line represents 
the luminosity expected on the assumption of the con- 
stant energy-conversion efficiency (rj). In model A, the 
luminosity is kept around the Eddington luminosity irre- 
spective of the mass accretion rate in the super-critical 
accretion regime, although the luminosity is in propor- 
tion to the mass accretion rate in the sub-critical accretion 
disk. Thus, the energy-conversion efficiency is remarkably 
reduced in the super-critical accretion due to significant 
photon-trapping effects, rj ~ rh' 1 . The calculated lumi- 
nosity in model B in the super-critical regime is larger than 
in model A owing to relatively ineffective photon trapping. 
Even so, the energy-conversion efficiency is also categor- 
ically on decrease in the super-critical accretion regime, 
rj cx rh' 1 / 2 . In both models, the results are in good agree- 
ment with the analytical estimates (§2). 

The dotted curve is the luminosity calculated based on 
the slim-disk model (Watarai, Mizuno, & Mineshige 2001). 
The energy is assumed to be mainly dissipated on the 
equatorial plane in the slim disks, however, the luminosity 
of the slim disk is not kept around the Eddington luminos- 
ity but increases with the mass accretion rate. Hence, the 
slim disks considerably overestimate the luminosity com- 
pared with our similar model (model A). The luminosity 
is also larger in the slim disk than in model B in the region 
of rh — 10 — 10 3 . The photon-trapping effects are already 
appreciable at L<L^, while they are only substantial at 
L >Le in the slim-disk treatment. The slim-disk model 

thus underestimates the photon trapping (see §1). 

Numerical results are roughly consistent with the ana- 
lytical assessments plotted based on equations (5) and (6), 
but deviations originate for the following reason. In nu- 
merical simulations, the energy transport at each height 
is actually solved by taking radiation transfer into consid- 
eration. On the other hand, the analytical assessments 
are merely obtained by the comparison of the accretion 
time-scale and the radiative diffusion time-scale. 

4.3. Effective temperature 

Since photon trapping tends to lower F s z urf at smaller 
radii, the disk with super-critical accretion rate has a flat- 
ter effective temperature profile in the vicinity of the in- 
ner edge. We calculate the effective temperature, T c g = 
[i^ urf /cr] 1//4 , and plot in Figure 4 the effective tempera- 
ture distributions for model A (solid curves) and model 
B (dashed curves), respectively. Although we adopt M — 
10M Q , the results can be applied to cases with other M, if 
we vary the absolute value of T e ff in proportion to M -1 / 4 . 
The effective temperature for low m, rh — 0.1 and 1, is 
proportion to r~ 3 / 4 as in the standard disk (Shakura & 
Sunyaev 1973; Pringle 1981), since the photon-trapping 
effects are not appreciable when the mass accretion rate 
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is sub-critical (see Figure 3). In the case of the super- 
critical accretion disk, the slope of the effective temper- 
ature profile becomes flatter within the trapping radius, 
though we still have T c g oc r~ 3 / 4 at r > r trap . For 
m > 10 3 in model A, especially, the profile within the 
trapping radius is almost flat and the maximum tempera- 
ture begins to decrease with farther rise in rh. This makes 
a marked difference from the prediction by the slim disk, 
which shows T c s oc r -1 / 2 in the high rh limit (Watarai 
& Fukue 1999; Watarai et al. 2000; discussed later) 
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Fig. 4. — The effective temperature profile for model A (solid 
curves) and for model B (dashed curves) for various values of m. In 
the case of fn < 1, the effective temperature is proportion to r~ 3 / 4 , 
as in the standard disk. However, the profile drastically changes in 
the super-critical accretion flow (m 2> 1), and becomes significantly 
flatter at the region of r < rtrap- In model A, the maximum tem- 
perature decreases with increase in the mass accretion rate and the 
profile is nearly flat when m > 10 3 . 

4.4. Spectral index 

As we have seen above, the effective temperature pro- 
file is very sensitive to the functional form of q v - ls {z). The 
vertical distribution of viscous heating rates would be ob- 
servationally distinguished when the luminosity is com- 
parable to or more than the Eddington luminosity. The 
spectral energy distribution (SED) is strongly affected by 
the effective temperature profile, which in turn depends 
on q v i s (z). We have calculated expected SED based on 
a superposition of blackbody spectra, B v \T c s{r)\, emitted 
from the disk surface, and plot in Figure 5 the spectral 
index, Q = dlogL„/dlogp, in the range between 0.1 keV 
and 1.0 keV, as a function of the luminosity normalized by 
the Eddington luminosity. Here, the circles and squares 
indicate the spectral indices for models A and B, respec- 
tively. It is known that the spectrum of the standard type 
disk but with T c g oc r~ p is composed of three compo- 
nents: Rayleigh- Jeans part (L v oc v 2 ) at low v, Wien part 
[oc i/ 3 cxp(— hv/kTcff)] at high v, and a power-law part 
(oc v z ~ 2 / p ) in the intermediate frequencies (Kato, Fukue, 
& Mineshige 1998), where v is the photon frequency, h 
is the Planck constant, and k is the Boltzmann constant. 



Thus, difference in p should manifest in this middle part. 
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Fig. 5. — The spectral indices as functions of normalized luminos- 
ity, -L/I/E- Here, we define the spectral index as £ = d log L„/ci log v 
at 0.1-1.0 keV. The circles and squares correspond to models A and 
B, respectively. The spectral index is sensitive to the distributions 
of the viscous heating rates in the vertical direction, thus it can be 
determined observationally. 

The effective temperature for rh = 0.1 and 1 is propor- 
tional to r -3 / 4 , however, the index deviates from £ = 1/3 
for p — —3/4 expected from the standard-disk model, since 
we see a Wien cut-off on the high-frequency side. The 
effective temperature profile in model A becomes flatter 
and the maximum temperature begins to decrease when 
rh > 10 2 , so that the peak of the SED shifts to the low- 
frequency side. Thus, the index drastically decreases in 
accordance with increasing rh, although the luminosity 
stays around the Eddington luminosity. In model B, on 
the other hand, the maximum temperature continues to 
increase with rise in the mass accretion rate. Therefore, 
the Wien cut-off is out of the range in the super-critical 
accretion. As shown in Figure 4, the effective tempera- 
ture profile in model B is flatter than in the standard disk, 
oc r~ 3 / 4 , but steeper than in the slim disk of high rh limit, 
oc r~ 1//2 . Therefore, the index is settled down between 1/3 
(standard disk) and —1 (slim disk), and the luminosity in- 
creases in accordance with rise in the mass accretion rate. 
We can thus obtain information of q vis (z) by the X-ray 
observations of luminous X-ray objects. 

4.5. Most luminous radius 

The radius, at which differential luminosity [(dL/dr)Ar] 
is maximum, is at (75/16)r g in the standard disk. This 
'most luminous radius' does not necessarily coincide with 
the inner edge radius, if the photon trapping is effective, 
since it is more prominent at small radii. We show the 
most luminous radii, r max , normalized by r g in Figure 6. 
The circles and squares indicate r max /r g for model A and 
model B, respectively. This figure clearly demonstrates 
for model A that the larger rh is, the larger becomes the 
most luminous radius, as long as ra > 10 2 . In model 
B, the most luminous radius stays at around the inner 
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edge of the disk, but the photon trapping significantly re- 
duces the overall radiative flux within the trapping radius. 
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Fig. 6. — The most luminous radii, at which (dL/dr)Ar is max- 
imum, normalized by r g are shown as functions of normalized mass 
accretion rate, m. The circles and squares correspond to models A 
and B, respectively. Since the photon-trapping effects are most ap- 
preciable in the vicinity of the inner edge, the most luminous radius 
shifts from the inner edge to the outer region at high fn in model A. 
In model B, in contrast, the innermost ring is still most luminous, 
although the radiative flux at the surface is reduced within rtrap 
(see Figure 3). 

4.6. Rapid accretion 

So far, since we adopt relatively large absorption coeffi- 
cient, the gas is approximately in equilibrium with radia- 
tion, T gas ~ T ra d, where T ra d is the radiation temperature. 
Thus, the energy dissipated by viscosity can be promptly 
converted into radiation energy, and that energy is trans- 
ported toward the disk surface by radiation transfer. How- 
ever, if the absorption coefficient is much smaller, as in the 
case that we adopt free-free absorption, it is possible that 
the energy dissipated by viscosity is kept as thermal en- 
ergy and is not converted into radiation energy. Therefore, 
the resultant luminosity must further decrease because of 
this effect in addition to the photon-trapping effects, if the 
gas accretes onto the black hole without achieving radia- 
tive equilibrium. This occurs, when a is relatively large 
a >0.03 (Bcloborodov 1998). 

To demonstrate this phenomenon, we employ the Rosse- 
land mean of the free- free absorption coefficient, Kg = 
1.7 x 10- 25 T g aI /2 (p/m p ) 2 (Rybicki & Lightman 1979), and 
calculate and plot the luminosity against the radial veloc- 
ity for model A (circles) and model B (squares) in Figure 
7, where the mass accretion rate is fixed to be m = 10 2 . 
It is found that the accretion disk gets much fainter at 
the regime of £ = v r /vg > 10~ 2 [corresponding to a>0.04 

if H/r ~ 0.5, since v r ~ a{H /r^v^}, where vg is the 
free-fall velocity. The accretion disks with huge accre- 
tion rate and large radial velocity tend to have extremely 
small energy-conversion efficiency and would be identi- 
fied as faint objects. Note, however, that this criterion 



of the radial velocity might vary, depending on the verti- 
cal density profile, since the free-free absorption coefficient 
strongly depends on the gas density, as well as the gas tem- 
perature. If the density and temperature of the gas is high 
around the equatorial plane, the thermal energy would be 
converted into the radiation energy there, and the disk 
might be somewhat brighter. Anyhow, the disk tends 
to be fainter in more rapid accretion in cooperation with 
the photon-trapping effects, although the detail examina- 
tion requires multi-dimensional radiation-hydrodynamical 
simulation. (This will be performed in future.) 
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Fig. 7. — The luminosity of the flow as a function of the accretion 
velocity, where the circles and squares indicate model A and model 
B, respectively. Here, we adopt the Rosseland mean of the free-free 
absorption coefficient. The accretion disk becomes dimmer in the 
regime of v r > W- 2 v H (£ > 1CT 2 ), since the energy dissipated by 
viscosity is not sufficiently converted into radiation energy but kept 
as thermal energy of the gas. The luminosity of the accretion disk 
is suppressed by this effect as well as the photon-trapping effects in 
the case of the disk with rapid accretion (i.e., large a > 0.04). 

5. DISCUSSION 

5.1. Comparison with the slim-disk model 

As a model for describing super-critical disk accretion 
flow, slim disks arc often utilized. However, we wish 
to note that the slim disk does not accurately treat the 
photon-trapping effects in the disk accretion case. The 
slim disk can describe the dynamical structure of the 
flow reasonably well, since the energy balance is given by 
Qvis ~ Qadv and does not depend on the radiative loss. 
The slim disk is not able to correctly estimate the radiative 
flux at the disk surface, in the sense that it overestimates 
the luminosity in all ranges in comparison with model A, 
and in the range of m between 10 and 10 3 in model B. 
Moreover, we found that the effective temperature pro- 
file becomes significantly flatter in model A, but slightly 
steeper in model B, T c g oc r~ p with p > 1/2. Hence, 
the profile strongly depends on the distribution of viscous 
heating rates. 

The deviations of our results from the slim disk orig- 
inate for the following reason. In the slim disk, using 



8 



Ohsuga et al. 



the vertically-integrated approximation, the accreting ve- 
locity of the flow is investigated by solving the equation 
of motion in the r-direction, and the radiative flux from 
the disk surface is determined as a function of radius 
so as to satisfy the vertically integrated energy equation, 
Q v i s = crT^ e + Qadv, where Q a dv is the advective cooling 
term. The resultant luminosity increases in accordance 
with rise in the mass accretion rate, and the effective tem- 
perature profile is obtained as T e ff oc r~ x / 2 in the limit of 
large rh. However, the radiative flux from the disk surface 
is assumed to be related to the temperature on the equa- 
torial plane as <tT^ s ~ oT c 4 /t. Since the radiation energy 
produced near the equatorial plane can not be radiated 
away immediately at the same radius but is advected in- 
ward at r < r trap , the effective temperature in reality does 
not always reflect the temperature of the equatorial plane 
at the same radius. The relation of oT e 4 ff ~ aT^/r actually 
holds, only when the photon-trapping is not appreciable. 
Hence, the slim disk cannot give the accurate radiative flux 
at the disk surface, although a part of the photon-trapping 
effects, the advection of the energy, is taken into consid- 
eration in the equation of energy balance. In the present 
study, on the other hand, we have actually solved radia- 
tive transfer in the z-direction (in the comoving frame) to 
investigate the radiative flux from the disk surface. Con- 
sequently, the observed luminosity and the effective tem- 
perature profile deviate from the prediction of the slim 
disk. 

5.2. Future work 

We have so far assume that the disk size is practically in- 
finite, but the super-critical accretion disk would be much 
dimmer when the disk size, r out , is smaller than the trap- 
ping radius. In model A, most of energy arises near equa- 
torial plane of the disk, and thus cannot arrive at the disk 
surface if r trap > r out . As a result, the luminosity in such 
a case must be extremely small, L <C Le- In model B, the 
luminosity is estimated as 



/Tout 

2 / 27rrQ v 

J fin 



= ^ Merit C 2 



2m \ 
3ft" J 
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(23) 



when r trap > r out . For example, this equation indicates 
L/Le ~ 3.1 and 0.11 for r out = 10 2 r g and 10r g , re- 
spectively, in comparison with L/L-g ~ 30 for the case 
of r out S> r trap , where we assume r- m = 3r g , rh = 10 4 , and 
h = 0.5. Such a situation might be realized in the accre- 
tion disks in the compact binaries, in the central regions 
of Type II supernovae (e.g., Mineshige et al. 1997) and in 
gamma-ray bursts (e.g., Narayan, Piran, & Kumar 2001). 

Throughout the present study, we describe the simple 
model for accretion motion without treating the hydrody- 
namical simulations. If the radiative flux force at the disk 
surface is stronger than the gravity of the black hole, the 
gas will be blown out like a wind (Eggum, Coroniti, & 
Katz 1987, 1988; Okuda, Fujita, & Sakashita 1997; Fujita 
& Okuda 1998). Since the ratio of the radiative flux force 



to the gravity in the z-direction is roughly estimated as 
~ m(r/r g ) _1 F s z urf /(5 v i s , the photon-trapping effects tend 
to prevent outflow by suppression of F s ^ rf /Q v is- Indeed, 
the gas would accrete onto the black hole without being 
blown out in model A (see Figure 3). Further study de- 
mands multi-dimensional radiation-hydrodynamical simu- 
lation. 

Moreover, we notice two very important improvement 
to be made. One is to include convection effect, which 
seems to occur when disk is radiation pressure-dominated 
(Shakura, Sunyaev, & Zilitinkevich 1978; Agol et al. 
2001). Even if the viscous heating is effective in the vicin- 
ity of the equatorial plane as in model A, the convection 
might more efficiently transfer the energy from the equato- 
rial plane to the disk surface than radiation. Then, resul- 
tant properties of the disk would become close to model 
B, but it is necessary for investigating in detail to solve 
two dimensional hydrodynamical equations coupled with 
radiation transfer. Second is to take the relativistic effect 
into consideration. It would play a significant role around 
the inner edge of the disk. 

5.3. Observational implications 

There seem to be several possible sites of super-critical 
accretion flow, where the photon-trapping effects may play 
a significant role. Firstly, we mention micro-quasars which 
compose a peculiar subclass of X-ray transients seemingly 
containing black holes (see Tanaka & Shibazaki 1996 for a 
review) . There is an indication that they emit rather high 
luminosity, near Eddington. In the case of GRS 1915+105, 
for example, its black hole mass determination was done 
as 18.6 ± 2.2M by Borozdin et al. (1999) based on the 
model fitting of the X-ray spectra. Independently, Greiner, 
Cuby, & McCaughrean (2001) obtained 14±4Af Q through 
observations of Doppler shifted CO lines of the secondary 
star. Thus, its X-ray luminosity of several times 10 39 erg 
s _1 at the peak is close to the Eddington luminosity of 
(2.3 ± 0.3) x 10 39 erg s _1 . Bursting behavior could be 
explained by relaxation oscillation between the slim and 
standard disk (Honma, Kato, & Matsumoto 1991; Ya- 
maoka, Ueda, & Inoue 2001). If this is the case, this pro- 
vides strong evidence for the presence of a super-critical 
accretion flow. 

Second one is ultra luminous X-ray sources (ULXs). 
ULXs with intermediate-massive black holes are succes- 
sively discovered in nearby galaxies and some of them seem 
to emit large luminosity, corresponding to the Eddington 
luminosity (Okada et al. 1998; Colbert & Mushotzky 1999; 
Mizuno et al. 1999; Makishima et al. 2000). In addition, 
an intermediate massive black hole of >10 3 M Q which is 

discovered in M82 (Matsumoto et al. 2001) is also a good 
candidate of super-critical accretion, although it is not al- 
ways clear if their luminosities are really close to the Ed- 
dington luminosities, since there is no accurate mass esti- 
mations. 

Narrow-line Seyfert 1 galaxies (NLSls), which are be- 
lieved to have a central black hole, are also considered to 
be sites of super-critical accretion (Mineshige et al. 2000). 
Since the NLSls exhibit large soft X-ray excess (Boiler, 
Brandt, & Fink 1996; Otani, Kii, & Miya 1996; Leighly 
1999; Hayashida 2000), mass of the central black hole is 
estimated asM~ 10 5_7 M Q by a comparison with soft X- 
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ray bump in galactic black hole candidates. Such moderate 
mass is also supported by reverberation mapping (Laor et 
al. 1997) and X-ray variability method (Hayashida 2000), 
and is consistent with narrow Balmcr lines of NLSls. 
Thus, the observed luminosity, L ~ 10 43 ~ 44 , is near the 
Eddington luminosity, and it implies that the accretion 
flow is super-critical. 

The final candidate is gamma-ray bursts (GRBs) (Pi- 
ran 1999, Meszaros, Rees, & Wijers 1999; Meszaros 2001). 
Interestingly, many models (such as mergers of double neu- 
tron stars, neutron star-black hole mergers, collapsars, 
and so on) predict a similar configuration as a final prod- 
uct; namely, a stellar-mass black hole surrounded by a 
massive torus with a mass of 0.01 - a few M Q . Then, 
gas accretion onto a new-born black hole may be the ori- 
gin of huge energy release and high Lorentz factor (e.g. 
Narayan, Paczynski, & Piran 1992; Narayan, Piran, & Ku- 
mar 2001). A simple estimation gives an enormous mass 
accretion rate. If gas with mass of M gas falls onto a black 
hole with a mass M on a time-scale of i acc sec, we find 
M- 2 x 10 33 (M gas /M Q )(i acc /ls)- 1 gs- 1 , that is, 

(24) 

(see, e.g., Narayan, Piran, & Kumar 2001). Then, the 
photon-trapping radius should be huge, r trap <~ 10 15 r g . 
Since the disk size is, at most, about the size of giant 
stars, r out - 10 13 " 14 cm - 10 7 ~ 8 (M/SMq)" 1 ^, it is very 
likely that the emission from the entire disk surface is to- 
tally blocked as we already mentioned, although neutrino 
energy loss may be essential in such a massive disk (Ruffcrt 
& Janka 1999; Popham, Woosley, & Fryer 1999). 

6. CONCLUSIONS 

By employing a simple accretion flow model, we have in- 
vestigated the photon-trapping effects in the super-critical 
black hole accretion flows by solving radiation transfer and 
energy equations of gas as well as radiation. The present 
results are summarized as follows. 

(1) The radiative flux at the disk surface is reduced 
within ~ mr g due to photon trapping. Then, the luminos- 
ity is kept around the Eddington luminosity if the viscous 
heating is effective only in the vicinity of the equatorial 



plane. Thus, the energy-conversion efficiency is drasti- 
cally on decrease, r\ = L/Mc 2 ~ to -1 . Even if the energy 
is dissipated equally in the disk, the energy-conversion ef- 
ficiency is also suppressed in accordance with rise in the 
mass accretion rate, r/ oc rnT 1 / 2 . 

(2) The slim disk can not accurately describe the effec- 
tive temperature profile, since it is assumed to be related to 
the temperature on the equatorial plane as T c ff ~ Tc/r 1 / 4 . 
As a result, the slim disk overestimates the luminosity of 
the disk at L <~ Le- The photon-trapping effects are al- 
ready appreciable at the luminosities as low as <Le, while 
they only become critical at L>3L^ according to the slim 
disk calculations. 

(3) Through the photon-trapping effects, the effective 
temperature profile becomes flatter than T e g- oc r~ 3 / 4 
within the photon trapping radius in the super-critical ac- 
cretion. If the energy is dissipated in the vicinity of the 
equatorial plane, the profile is almost flat (T c ff oc r~ p with 
p ~ 0) and the maximum temperature decreases in accor- 
dance with the rise in the mass accretion rate. 

(4) The information regarding the distribution of the 
viscous heating rates, q v - ls (z), can be obtained through the 
observations of luminous X-ray objects at L>L^- In ac- 
cordance with rise in the mass accretion rate, the spec- 
tral index, C = d\ogL v /d\ogv [0.1 keV - 1.0 keV] for 
M = 1OM , increases but in the different ways, depending 
the functional form of q v i s (z). When the gas is heated up 
uniformly, independently of the vertical height, especially, 
the spectral index stays between —0.5 and 0, although the 
index must be further small, <C — 1, in the case that the 
energy is dissipated around the equatorial plane. More- 
over, the most luminous radius shifts to the outer region 
as long as to ;§> 10 2 in the latter case. 
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